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ABSTRACT 

Photoionization of the Ce'^+ - Ce^+ process has been studied using the ran- 
dom phase approximation with exchange method in the energy region 100-150 
eV. Comparison of our results with the recent measurement [Miiller et al, Phys. 
Rev. Lett. 101, 133001 (2008)] confirms the suppression effect of the carbon 
cage in the endohedral fuUerene Ce@C^2 photoionization. The reasons for the 
cause of the confinement resonance and the suppression effect are discussed. 



1. Introduction 

The motivation to study the photoionization of the Ce^+ - Ce*+ process is 
the disturbing discrepancy between the theoretical calculations and the experi- 
mental measurements of the photoionization of an endohedrally confined atom. 
The theoretical calculations [1-6] indicate the presence of strong confinement 
resonances for the endohedral fuUerenes, such as XeQCeo [6]. However, the 
experimental results [7-11] demonstrate a great suppression of the photoioniza- 
tion cross section of an atom encapsulated in the carbon sphere, for example 
the photoionization of Ce@C82 [7], Pr@C82 [8], Ce@C^2 [9> 10] > and Dy@C82 
[11]. 

Both experiments [12, 13] and theoretical study [14] have indicated that 
the Ce atom in the endohedral fullerene Ce@C82 is located at an off-centered 
position adjacent to the carbon cage and the encapsulated Ce atom donates 
three valence electrons to the carbon sphere. The electronic state of Ce@C82 can 
be formally described as Ce^+QCg^ [13]. The same charge state can be found for 
the Ce atom in the endohedral molecule Ce@Cg'2 [10]. Since the photoelectron 
ionized from the Ce^+ ion of the Ce@C82 molecule will be multiply reflected to 
different directions by the carbon cage, it is difficult from a theoretical point 
of view to consider the photoionization of an off-centered Ce^+ ion. However, 
the photoionization of the Ce^+ ion can be calculated by using our recently 
developed random phase approximation with exchange (RPAE) method [15] if 
the necessary modifications in the computer code are made. The results should 
further confirm the cage suppression effect if agreement is obtained with the 
experimental data for the Ce^+ - Ce^+ photoionization [10] or a transparent cage 
model can be set up for the photoelectron if it agrees with the photoionization 
cross section of the Ce@Cg2. 

The RPAE method, which allows for the inclusion of both intrashell and 
intershell correlations has been developed recently by Chen and Msezane for 
atoms(ions) with an outer open-shell [15], or with an inner open shell [16] and 
has been successfully used to study the M — ef photoionization of the ions 
Xe"*" [15] and 1+ [16] and the photoelectron angular distribution asymmetry 
parameter (3 of the Sc 4s electron [17]. In this calculation our RPAE codes are 
modified to include the intershell coupling between the Ce^"*" 5s — ep, 5p — es, d 
and Ad — Af transitions. A new computer code is used to study the Ce^"*" 
photoionization. 

2. Theory 



The RPAE equation and the symbols and operators in the equation for an 
atom with an outer open-shell is given by Eq. (1) of Ref [15]. Similar terms 
will be added for the switch of other electron pairs. All the matrix elements 
have been derived and presented in the Appendices of Refs. [15] and [17]. The 
Coulomb matrix elements which are needed to evaluate the intershell coupling 
between the 5s — ep, 5p — es, ed and 4d — 4/ transitions of the Ce^+ ion, are 
given in the Appendix of this paper. 

The ground state of the Ce^"*" ion has the configuration [Xe]4f(^F). Here we 
assume the 5d and 6s^ electrons of the Ce atom, which are located outside the 
4/ orbital have been transferred to the carbon cage in the endohedral fullerene 
Ce@C^2- 

Since we try to include the intershell couplings among the 5p — es, ed ,5s — 
ep, and 4/ — ed, eg transitions , the following combined core with discrete and 
continuum electron states have been included in the calculation. 

We have a total of 21 channels from the 5p + hu ^ ed, es transitions to the 
states: 

M^%s^5pHf{^D, iG, 3£), 3G, ^F,)edCD, ^F, ^G), 
M^'^5s^5pHfCD, ^D)es{^D), 
M^°5s^5pHf{^F, ^F)es{^F), 
4^105525/4/(10, ^G)esCG), 

a total of 6 channels from the 5.s + hv ^ ep transition to the states: 
4^105525/4/(1^, ^F)epCD, ^F, 2G), 

and a total of 2 channels from the Af + hu ^ ed, eg transitions to the states: 

M^°5s^5p'^{\S)edCD), 
Ad^'^5s^5p'^CS)egCG). 

The closed shells of Is, 2s, 2p, 3s, 3p, 4s, and 4p axe not listed above. As the 
photoionization process is caused mainly by the autoionization of the 4/ subshell 
our code has several special subroutines to treat the M^^Af + hu ^ Ad^if^ 
transition, with a total of 14 channels: 

4^105^25^64/ + /jj, ^ 4d^5s''5p%4fCH))CF, 2G) (1) 

4d^°5s'^5p^4f + hu ^ 4dF5s'^5p^{4f{^F))CD, ^G) (2) 

Ad^^hs^bp^Af + hv ^ Ad^hs^bp^{Af{^P)){^D, ^F) (3) 

4di°5s25p^4/ + hy^ 4d^5s''5p%4fCl))CG) (4) 

4d^°5s'^5pHf + hu ^ 4d^5s'^5p'^{4fCG, ^D))(^D, ^F, ^G) (5) 



The Ce^+ ground state and the core wave functions were obtained through 
the self-consistent Hartree-Fock (HF) calculation. Then the radial functions of 
the discrete and continuum electron were obtained by solving the linear HF 
equations without self-consistency using those core wave functions. Each radial 
function has been represented by 2000 points. After evaluating the dipole matrix 
elements and the Coulomb matrix elements of the time-forward type and the 
time-backward type, the RPAE equation was solved to obtain the partial cross 
sections with a total of 15 states, 15 states, and 13 states. All types 
of matrix elements are evaluated using the equations found in the Appendices 
of Ref. [15] and Ref. [16]. Equations for the Coulomb matrix elements of the 
intershell coupling between the Ce^+ 5p— es, ed, 5s — ep and 4d — 4/ transitions 
are found in the Appendix of this paper. 
Results 

Its well known that there are two kinds of giant resonances [18]. One is the 
shape resonance in which the properties of the giant resonance is determined 
mainly by the effective potential for the / electron. The overlaps of bound 
/ orbital with the 4rf orbital are very small. Most of the 4d — > / oscillator 
strength is associated with the continuum state. The photoionization of the 
4d — ef channel in the atoms Xe, Ba, and I [19] and ions Xe+ [15] and 1+ 
[16] are the examples of this type of giant resonance. The resonance results 
from a one step process: continuum enhancement due to the centrifugal-barrier 
shape resonance. The second type of giant resonance corresponds to a decaying 
discrete resonance, which results from a two-step process: photoexcitation of 
a 4d electron into the 4/ subshell, followed by autoionization of the 4d^4/^+^ 
state where N is the initial occupation number for the 4/ electron. An example 
of this type of giant resonance is found in the photoionization of the rare earths 
elements. The decaying discrete resonance is related to the so-called collapsed 
/-wave function. In this situation the inner well is deep enough to support a 
bound state and the 4/ orbital "collapses" into the inner well. 

The photoionization of the Ce^+ id electron belongs to the second type of 
giant resonance. The mean radus of the 4/ electrons is only 0.96 a.u., which 
is much smaller than that of the 5s electrons (1.58 a.u.) and 5p electrons 
(1.75 a.u.). Most of the Ad ^ f oscillator strength is associated with the 4/ 
state. We performed two calculations: one with the processes represented by 
equations (l)-(6), and the other without these processes. The results show 



that the photoionization cross section without the Ad — 4/ transition is about 
only 10% of the results when these transitions are included. It demonstrates 
that the photoionization cross section is mainly (90%) caused by two processes. 
The first process is to photoexcite the id electron to the 4/ subshell, then the 
autoionization of the 4/ subshell causes the photoionization of the 4/, 5s, and 
5p subshells. 

Figures 1 gives our calculated photoionization cross sections for the Ce^+- 
Ce^+ transition in the energy region 100-150 eV. Dotted, dashed and dash- 
dotted curves represent, respectively the partial cross sections for the ^F, 
and states. The first peak, located at 101.7 eV, is contributed mainly (88%) 
by the partial cross section of the '^F state. The most important photoexcited 
states arc 

4d95s25/(4/2(iG))(2F), and 

4^95525/(4/2 (3^))(2^)_ 

The most important final states through intershell coupling among the 5p—ed 
5s — ep and M — 4/ transitions are 
4di05s25p54/(iG)ed(2F), 
Ad^%s^bp^Af{^D)edCF), 
4d^%s^5pHfi''G)edCF), 
4d^°5s5p^Hf(:^F)ep(^F), 
4d^°5s5pHf{^F)ep(^F). 

Other transitions will also contribute to the cross section but to a lesser 
extent. We note that the above five states contribute 86.0% to the partial cross 
section of the ^F state in the first peak at 101.7 eV. 

The second peak is at 118.0 eV with a maximum cross section of 47.9 Mb. 
The 2£), '^F and states contribute respectively 2.4 Mb, 3.9 Mb and 41.6 Mb 
to this peak. The transition to the state accounts for 86.8% of the total cross 
section. Therefore, the most important processes are first Ad electron excited to 
the states 

4d95s25/(4/2(3H))(2G), and 

4d95s25p6(4/2(lj))(2(5)_ 

Then autoionization follows to the final states 

4di05s25p6(i^)eg(2G), 
4di05s25p54/(3D)erf(2G), 
4di"5.s25/4/(3F)erf(2G), 
4d^°5s5p^Af{^F)ep(^G). 



Since the 4/ — eg transition has the cross section 20.4 Mb, which is almost 
half the cross section (41.6 Mb) of the state, this transition is the most 
important process for the second peak. 

The third peak has a maximum of 57.0 Mb and is due to the partial cross 
sections for the "^D, and states whose values are respectively, 23.7 Mb, 
26.0 Mb, and 7.3 Mb. Unlike in the first and second peaks, for the third peak 
the and states contribute significantly, while the state contributes 
moderately to the cross section. The most important final states are 

Ad'^^5s^5pHf{^G)ed{^D), 

4d'^°5s^5p^4f{^F)ed{^F), 
M^°5s^5pHf{^G)edCF) 
4di°5s5/4/(3F)gp(2F), 
4di05s25/(i5)e^(2G). 

From the above analysis we can see that the ^D, ^F and '^G terms of the 
core wave function are more important than the ^D, ^F and ^G terms. 

Figure 2 displays the comparison of our RPAE results with the experimental 
data. Solid and dotted curves represent, respectively our calculations and the 
data from the measurement [10]. The agreement is reasonable. The two main 
peaks in the experiment are closer to each other than those in our calculation. 
Both the theoretical and experimental results are much larger than the cross 
section for photoionization of Ce@Cg2 (see section a of Figure 2 in Ref. [10]). 
The photoionization cross section for Ce@Cg'2 has a maximum of 20 Mb around 
123 eV. This value is much smaller than that of our third peak 57.0 Mb at 123.2 
eV. The reduction eflfect from the carbon cage can also be seen from comparison 
of our results with the experimental data of the endohedral fullerene Ce@C82 
[7]. The photoionization cross section of the 4d — 4/ giant dipole resonance 
in Ce@C82 is estimated to be 14.3 Mb at 130 eV, which is also smaller than 
our results of 20.9 Mb at 130.0 eV. Therefore our calculation further confirms 
the suppression eS^ect obtained by the experiment when an atom encapsulated 
inside the carbon cage is photoionized. 

Many theoretical calculations demonstrate confinement resonances [1-6]. How- 
ever, all the experiments show no confinement resonances. The measurements 
demonstrate a great suppression of the photoionization cross section by the car- 
bon cage. It is hypothesized in Ref. [10] that additional decay channels for the 
Ce Ad vacancy may exist. After analyzing carefully the theoretical models and 



the experimental arrangement, we are led to the conclusion that the discrepancy 
between the theory and measurement may be due to the assumed location of the 
atom inside the Ceo by theory. The location of the atom in the calculation of 
the confinement resonances is always assumed to be at the center of the carbon 
cage. The photoelectron ionized from the atom will be reflected by the carbon 
sphere. For the atom assumed to be at the center of the Ceo the reflected wave 
and the incoming wave easily interfere with each other and create the resonance. 
However, the atoms in the current experiments involving the photoionization of 
CeQCgj, Ce@C82, Pr@C82 etc. are located at off-center positions and adjacent 
to the carbon cage. This geometrical configuration causes the photoelectron 
ionized from for example, the Ce^+ ion to be multiply reflected to different 
directions by the carbon cage. The reflected wave and incoming wave do not 
necessarily readily interfere with each other, but are easily absorbed by the car- 
bon cage in the process of photoionization. This may be the reason for the great 
suppression observed in the experiment. 

To resolve directly the long-standing discrepancy between the measurements 
on the one hand and theoretical predictions on the other on the photoionization 
of an endohedral fuUerene, we recommend that both measurements and calcu- 
lations be performed on the photoionization of A@Ceo (atom at the center of 
Ceo cage). 
Conclusion 

In conclusion, we have performed a RPAE calculation for the photoioniza- 
tion of the Ce^+ - Ce^"*" process. The reasonable agreement in magnitude and 
shape with the recent measurement [10] confirms the suppression effect of the 
carbon cage in a endohedral fuUerene. The plausible reasons for the confinement 
resonance or suppression effect have been advanced and discussed as well. 
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Appendix: Coulomb matrix elements for the intershell coupling between Ce^"*" 
5s — ep, 5p — es, d and Ad — 4/ transitions 



In the following equations Z"^ is an open shell and the other Z" subshells are 
closed shells. G^J^J etc. is the fractional parentage coefficient. 
(1) For the Coulomb interaction 

\iT-\L'^s'^iT[L2S2]iT^\L\s[]rs^^ >^ 

\ir [LsSs]ir-'[L'^Si,K^ [L,S,] [L'XlkL'S' > 
the Coulomb matrix element of time forward type is: 

^L[s[ [L',,S'„L[,S[]'/^L']V^^r+T \ h h k 




1/2 Si S'^ 
1/2 S' Si 





The exchange part of the Coulomb matrix element of time backward type 



is: 



* < 



h h k 

* < kWC'Wh >< hWC'Wh > i?fe(ZiZ6;Z2Z3)(-i)'^""+'+"^+'^"+'^+'^-'' 



5*^ 1/2 S' 
5c" 1/2 S 



S[ S"' 1/2 
Sc" Si 1/2 



Lc" Li li 



where Lc"5c" are from coupling of l^'^^ii^S^]/;"^ [Li5i] [i^S*^]. The symbols () 
and { } are the so-called 3j symbol and 6j symbol, respectively. The symbol 
{ } with three lines and three columns is the so-called 9j symbol. C'^ is the 
normalized spherical harmonics, 



h k I3 

L' a L 



and 
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Figure Captions 

Fig. 1. The cross sections for the Cc'^+-Ce^+ photoionization process in 
the energy region 100-150 eV calculated in the RPAE approximation. Dotted, 
dashed and dash-dotted curves represent, respectively the partial cross sections 
for the symmetries ^_F. and ^G. The solid curve is the total photoionization 
cross section, which equals the sum of these partial cross sections. 

Fig. 2. Comparison of our RPAE results with the experimental data for 
the Cc''+ - Ce''+ photoionization process [10]. Solid and dotted curves are 
respectively, the calculated and measured results. 
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